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The precise knowledge of the beam energy, -Ebeam, at the LEP collider is important to reduce 
the systematic uncertainty on the W mass. The measurements by the LEP energy group 
can be cross checked using Z recoil events. Preliminary results of the four LEP experiments 
ALEPH, DELPHI, L3, and OPAL are presented. The combination of the results shows no 
significant deviation of the i?b cam value obtained using Z recoil events from the measurement 
by the LEP energy group. 



1 Determination of the LEP beam energy 



The very precise method to determine the beam energy that was used for the Z run of LEP with 
an uncertainty of about 1 MeV, the resonant depolarization method^, is not available for beam 
energies above 60 GeV because the electron and positron beams remain unpolarized beyond this 
energy. To determine the energy of beams of a higher energy - more than 100 GeV in the last year 
of LEP operation - another method had to be used. Nuclear magnetic resonance (NMR) probes 
were installed in 16 of the LEP dipole magnets. The beam energy is obtained using the magnetic 
field B measured by these probes: -Ebeam ~ § B ■ dl. This value has to be calibrated against the 
energy from the depolarization method at lower energies by assuming a linear relation between 
-f'beam and B to extrapolate to higher energies. The extrapolation method is cross checked using 
flux loop measurements, a spectrometer, and spin tune measurements^. The uncertainty on the 
LEP beam energy measurement from these methods is about 20 MeV. The LEP energy model 
also describes the time-dependent energy variation at each interaction point. An independent 
cross check can be obtained by measuring the kinematics of Z recoil events as described below. 
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Figure 1: The Feynman diagram for Z recoil events. The beam electron or positron emits a hard initial state 
photon. An on-shell Z boson is produced, which decays to a fermion-antifermion-pair. 



2 Z recoil events 



Figure ^ shows the Feynman diagram of the "radiative return" events used for the beam energy 
determination. Either the incoming electron or positron emits a high energy photon before 
the collision. The Z exchange process has a high cross section at energies around the Z mass. 
Therefore the energy of the photon is such that an on-shell Z boson is produced in most cases. 
In general, the photon escapes along the beam direction, but it is sometimes observed in the 
calorimeters. 

Due to photon emission, the effective centre-of-mass energy, v 7 ?, for fermion-antifermion-pair 
production is reduced from the nominal centre-of-mass energy \fs. The reduced centre-of-mass 
energy can be extracted from event kinematics using 
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Figure 2: Distributions of the effective centre-of-mass energy, vP, as measured by OPAL, for /i + fx~ (7) final 
states (upper left), r + r~(7) final states (upper right), and e + e~(7) final states (lower plots). Due to the large 
t-channel contribution, the recoil peak for the electron case is not seen as clearly as for the other channels. The 
"comb" structure at the end of the spectrum is due to the integral over different energies. 
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Figure 3: Effective mass of two-jet events, taken with the L3 detector, after kinematic fit, in the recoil peak 

region. 

where {9j) 6f is the angle between the (anti-)fermion and the photon. Alternatively, a kinematic 
fit can be applied to the event, imposing four-momentum-conservation. This fit significantly 
enhances the resolution of the spectrum of the effective mass m e fj. 

All four LEP experiments use the muon pair final state for this analysis. This channel has 
a very clear signature and low background. OPAL also uses the electron and the tau pair final 
states. Figure shows y / s'-distributions for lepton final states, as measured by OPAL. The 
radiative return peak at effective centre-of-mass energies of the Z boson mass is clearly visible. 
In addition, DELPHI, L3, and OPAL use two jet final states, making use of the high statistics 
of this channel. Figure El shows the distribution of the effective mass after kinematic fit for the 
L3 hadronic data. 

The analyses presented her d I I I I are performed on the 1997 to 2000 data sets with a total 
integrated luminosity of about 685 pb _1 per experiment. 

3 Beam energy extraction 

Different methods are used to extract the beam energy from the y/s' or m e a spectra. In general, 
the precise knowledge of the Z boson mass, mz = 91.1876 ± 0.0021 GeV"J, allows to measure 
the beam energy using the position of the radiative return peak. In the analyses either E^ eiim is 
assumed to be the nominal value, then a new mz value is determined, and from the deviation 
of the new Z mass from the precision measurement a shift in the beam energy can be extracted. 
Alternatively, the Z boson mass is assumed to be the precision mass, and the beam energy is 
extracted directly. 

For ALEPH muon pair events, a maximum likelihood fit is performed on the distributions 
of \fs! and on the "raw" effective mass 

m raw = 2 -VW2 ■ (1 - cos 012 ), (2) 

where pi are the measured momenta of the muons and 0i2 is the angle between them. 

From a likelihood fit to both spectra, Data - Monte Carlo differences AyV and Am raw are 
extracted. These values are used to create a \ 2 function, which depends on the beam energy. 
By minimizing \ 2 = /(-^beam)) a possible deviation 
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Figure 4: The "raw" effective mass (ieft) and the effective centre-of-mass energy (right) of muon pair events, as 
recorded with the ALEPH detector. Dots represent Data, whereas the histograms show Monte Cario predictions 
for different Z masses. The fines show the functions used in the maximum likelihood fit. 



from the nominal beam energy is obtained. 

In the second method, used by DELPHI and OPAL for lepton pair events, the Z boson mass 
is assumed to be the precision mass. The deviation from the nominal centre-of-mass energy, 
A-Ecms = 2 • AE'beam is expressed as 

A£ C MS = — - V^LEP, (4) 

x 

where function of the (anti-)fermion-photon angles only. A possible deviation from 

the nominal beam energy, AEbeamj is then extracted from comparing the results from a fit of a 
Breit-Wigner-like function to the Ai^cMS distributions for Data and Monte Carlo prediction as 
shown in Figure 03 

In an ansatz used by DELPHI for hadronic final states and by L3 for hadronic and muon 
pair events, an unbinned maximum likelihood fit, using Monte Carlo event weights 



/(vV,m test ) 
f{Vs',m z ) 



where /(vs',m) is a Breit-Wigner-like function and mt es t is the mass to be tested, is performed 
to the y/p or m e g spectra to extract the Z boson mass, ml*. Any deviation from the precision 
mass is attributed to a deviation of the beam energy from the nominal value using 

A^beam = "-^beam ~ ■ ( 6 ) 
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In the last method described here, used by OPAL for hadronic events, a Breit-Wigner-like 
function is fitted to the Monte Carlo y/s' spectrum to obtain the peak position m* mc . The same 
function is fitted to the data spectrum, leaving only the normalisation and the peak position, 
m Data' as function of .Ebeam, as free parameters. A distribution of the data peak position versus 
the beam energy is obtained, and the beam energy is extracted by requiring 



m Data(^beam) — ^mc- 



(7) 
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Figure 5: A_Ecm distributions from DELPHI, for Monte Carlo prediction (left) and /x-pair Data (right). A 

Breit-Wigner-like function is fitted to both spectra. 



4 Systematic uncertainties 

The main systematic uncertainties for lepton pair events are the determination of the lepton 
angles, the ISR modeling (studied by using different orders of photon radiation as implemented 
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in the KLfC2f^ Monte Carlo generator), and fit properties. Further uncertainty sources include 
background uncertainties, Monte Carlo statistics, and LEP parameters. The total systematic 
uncertainty for the beam energy shift observed in leptonic events is 16 MeV. 

For hadronic events, the systematic uncertainty is dominated by hadronisation uncertainties. 
They are studied by comparing the ARIADNE^, HERWIG^, and PYTHIaEI hadronisation 
models. Further major contributions are jet measurements and detector modeling. For the 
combined result from hadronic events, the total systematic uncertainty amounts to 43 MeV. 

For the combined result, the main contribution to the systematic uncertainty, amounting to 
15 MeV, is due to hadronisation. The total systematic uncertainty is 22 MeV. 



5 Results 

Tabled shows preliminary results of each experiment for leptonic and hadronic final states. A 
graphical representation is given in Figure H3 The final result, obtained by averaging the results 
from lepton pair and hadronic channels of all experiments, taking all correlations into account, 
yields 

A£ bcam = EgZ ~ E bfL = ("20 ± 25 ± 22)MeV. (8) 

This result demonstrates the good agreement between the beam energy measurement from 
radiative events and the standard LEP energy determination. 
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Figure 6: Graphical representation of preliminary -Bbeam measurements performed by the four LEP experiments. 
Shown is the shift from the nominal beam energy. The error bars indicate the total error. 



Experiment 


Channel 


A£ bcam [MeV] 


ALEPH 




-167 ± 91 ± 48 


DELPHI 


M + M~7 


+113 ± 75 ± 27 


DELPHI 


qq7 


- 55 + 53 ± 65 


L3 




+ 10+115 + 22 


L3 


qq7 


- 46 ± 33 ± 51 


OPAL 


e + e~7 


+ 40 ± 136 ± 78 


OPAL 




- 51 + 84 + 22 


OPAL 


r + r~7 


+301 ± 199 ± 148 


OPAL 


qq7 


- 66+ 34+ 70 


Combined 


l + l- 7 


+ 5+41+16 


Combined 


qq7 


- 52 + 24+ 43 


Combined 


ff 7 


- 20+ 25 + 22 



Table 1: A.Ebeam results from each experiment. The first error value represents the statistical contribution, the 

second the systematic contribution. 
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